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I try to present a small view of the properties and issues related to astronomical 
interferometry observations. I recall a bit of history of the technique, give some 
basic assessments to the principle of interferometry, and finally, describe physical 
processes and limitations that affect optical long baseline interferometry and which 
are, in general, very useful for everyday work. Therefore, this text is not intended 
to perform strong demonstrations and show accurate results, but rather to transmit 
the general "feeling" one needs to have to not be destabilised by the first contact 
to real world interferometry. 
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1 What is optical / IR long baseline interferometry ? 



1.1 In the beginning 



The d iscovery of the wave-property of light was probably made by lYoung 
( II8OOI ). who managed to produce light interfer e nces b y letting it go through 2 
close holes. Soon after this discovery, iFizeau I (Il85ll ) proposed and IStephan 
( 118741 ) tried unsuccessfully to use this wave-property of light to measure the 
apparent diameter of stars at the Observatoire de Marseille with a 80 cm 
telescope. 



Then, iMichelson and Pease I ( 1l92ll ) managed to measure the diameter of a 
star, Betelgeuse, equal to 0,047" with a relative accuracy of 10% using a larger 
interferometer (see Fig. 1, left) This experiment was very hard to carry out, the 
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first results (jPease I . Il921al : iMerrill I . Il922[ ) were so enco uraKing that a larger 
and more sensitive stellar interferometer (jPease I . Il93ll ) was built. However, 
technical drawbacks and mechanical instabilities were the major problems of 
these interferometers. Then, the project was abandonned during the second 
world war. 



When the very first radio-telescopes put into operation, (jReber &: Greenstein 



19471 ). the idea to coherentl y cornbine t ogether several antennas was used to 
build radio-interferometers (JRyle I . ll952l). This e nabled an increase in the spa- 
tial resolution of radio-collectors (jSmith I . Il952l ). An easy access to the mea- 
sures allowed a swift progress of th e instruments. This led to the VLA (Very 
Large Arrav. iButler et al. I . 120061 ) and to the VLBA (Very Long Baselines 



Array, I Cohen et al. I . Il975l ) 
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Fig. 1. Left: The interferometer of Michelson on the 100-inch telescope on the 
Mount W i lson ( taken from reference iMichelson and Pease I . Il92ll ). Right: The first 
Labevrie I ( 19751 ) interferometer (diagram is from the article). 



1.2 Given up and reborn ... 



In the optical domain, the techniques developed by Michelson were forgot- 
ten for more than 30 years . Duri ng this period, only intensity interferometry 
(jHanbury Brown fc Twiss I . ll956l ). which indirectly recombined the light, was 
developped and used for astronomy. 



The tr ue reborn of optical interferometry occurred in 1975 when A. iLabeyrie 
(I1975I ) managed to produce interference fringes on a star using two separate 
telescopes (see Fig. 1, right). 



This achievement was repeated in numerous institutes where many prototypes 



were built and then used for science (see Table 1). The initial c oncept was irn- 
proved with the use of delay lines and stationary telescopes (JMerkle I . Il986l ) 
rath er than moving optical tables a nd telescopes. The advent of spatial filter- 
ing (jCoude du Foresto et al. 1 . 119971 ) and sim ultaneou s photo metric calibration 
allowed one to get rid of systematic effects (jPerrin I . l2003al Jbl) . This enabled 
an advancement in accurate measurements. 



Table 1 
Working interferometers in the world in alphabetical order. NIR means Near Infra- 
Red (1-2. 5^m) and MIR means Mid Infra-Red (8-13/um). 



Name 


T. 
Diameter (m) 


slescopes 
# Combined 


Total 


Max. baseline 


A 


CHARA 


1 


2 


6 


330 


Visible, NIR 


COAST 


0.4 


3 


6 


47 


Visible, NIR 


ISl 


1.65 


3 


3 


85 


MIR 


Keck-1 


10 


2 


2 


85 


NIR, MIR 


MlRA-1 


0.25 


2 


2 


30 


Visible 


NPOl 


0.12 


6 


6 


64 


Visible 


PTl 


0.4 


3 


3 


110 


NIR 


SUSl 


0.14 


2 


2 


640 


Visible 


VLTI 


8/ 1.8 


3/3 


4/4 


130 / 200 


NIR, MIR 



1.3 Inter ferometry today 



Today's interferometry follows two different tracks: larger baselines to get 
higher angular resolutions and larger telescope diameters to reach higher mag- 
nitudes. The Keck-I and the VLTI, which combines both large baselines and 
large apertures, are leading the way. 



Keck-I: The Keck-I is a US project consisting of two very large 10 m seg- 
mented telescopes separated by 85 m. The Keck telescopes work in two main 
modes: independently or combined in the interferometric mode. This last mode 
uses two instruments: a coaxial re-combiner and a 2-telescope nuller. Today's 
K-band limiting magnitude is 10, and it can work in low (R ~ 20) and medium 



(R ^ 200) spectral resolutions. A differential mode and a phase referencing 
mode are foreseen for the future. However, the initially planned additional 
small telescopes for the interferometers (the "outriggers") have been aban- 
doned, limiting the Keck-I project to a 2 telescope only interferometer. 



VLTI: VLT (or Very Large Telescope) is European project mainly made of 
four 8 m telescopes (Unit Telescopes or UT) and a series of "small" 1.8 m tele- 
scopes (Auxiliary Telescopes or AT). Today, 4 ATs and 4 UTs are operational 
with baselines ranging from 16 m to 130 m (Fig. 2). Like the Keck telescope, 
the VLT can operate each 8 m telescope individually or combine up to three 
telescopes together using the VLTI or Very Large Telescope Interferometer. 3 
instruments are now operational at VLTI: 

VINCI: test instrument, 

MIDI: mid-infrared instrument using 2 telescopes 

AMBER: 3 telescope instrument in the near infrared. 

Today's VLTI is limited by vibrations, which are under investigation and 
should be solved in the next few years. The future of VLTI consists of a phase- 
referencing facility, PRIMA, and a series of projects for second-generation 
instruments (MATISSE, VSI, and GRAVITY). 



1.4 Definitions 

In this section, I will use the following notations: 

• t is time, 

• r is delay, 

• cj is the light pulsation and A is its wavelength: u = 2710/ \ 

• a; is the (x, y, z) coordinates vector, 

• "p^ is the (x, y) coordinates vector in the pupil plane, 

• s is the (x, y) coordinates vector in the detector plane, 

• M is the (m, v) coordinates vector in the Fourier plane. 

Other notations will be defined as they appear. 



Electromagnetic wave: An electromagnetic wave ( i? ( x , t), 5 ( x , t)) is a 
particular case of an electric and magnetic field that propagates in space with 
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Fig. 2. Aerial view of the Paranal mountain and the Very Large Telescopes (photo 
credit: Gerhard Hiidepohl) under-plotted on the VLTI stations. Small stations cor- 
respond to ATs and larger encircled ones correspond to UTs. 

time. The propagation of an electromagnetic wave is described by the Maxwell 
equations, and using the Lorentz gauge in space, it can take the following form: 



E{'x,t)=Eo{'^)e"^' 



;i) 



Light intensity: An optical or infrared detector is sensitive to the light 
intensity /(V); i.e. the time-average of the squared modulus of the electro- 
magnetic field at the measurement point: 



E\ 



■ t) 



(2) 



In this equation, s are the spatial coordinates in the detector plane (see 
Fig. 3), t is the time, and E is the electromagnetic field. 



star at infinity 



plane wavefront 



Pupil 



Screen 




S 




Fig. 3. The light propagation of an electromagnetic wave from a star at infinity 
to the detector plane (screen). Interferometry is interested in the study of the in- 
teraction between waves coming from two points of the initial wavefront. Notations 
come from the text. 



1.5 What does an interferometer measure ? 



Coherence between 2 waves: When observing an object with a given 
instrument, the light intensity / is the result of the superposition of many 
electromagnetic waves coming from different points of the instrument's pupil 
plane (see Fig. 3): 




(3) 



This expression applies when one is only interested in a 2-wave interaction: 



Ks 



E{^,t-Ti) + E{^,t-T2] 



(4) 



If one defines the mutual coherence function Fi 2 by: 



T^,2ir) = (E{^^,t-r)E*i^2,t) 



(5) 



The light intensity can be expressed in a simple way: 



/(^) = ri,i(o) + 12,2(0) + ri,2(r2 - n) + ri,2*(r2 - n] 



(6) 



Then one can define the complex coherence degree 71,2 (t) as: 



7i,2(r) 



ri,2(r) 



ri,i(o) + 12,2(0) 



(7) 



Light intensity becomes: 



/(^) = M^) + /2(^)] [1 + K (7i,2(r))] (8) 

which, in the case of wave-fronts at infinity, as described in equation 2, and 
considering the on-screen coordinate x and the wavelength A = 27rc/u, the 
light intensity becomes: 



J(^) = [/l(^)+/2(^)] 



2'ITX 
1 + /i COS ( — h 



(9) 



fj, being the modulus of 71,2(0) and (f) its phase. 



This co-sinusoidal modulation of the light is called an interferogram and is 
what one can measure using an interferometer [x is then modulated spatially 
- multiaxial instrument - or temporally - coaxial instrument - see Fig. 4). 





Fig. 4. Coaxial vs multiaxial beam recombination. In co-axial recombination, the 
fringes are scanned by a temporal modulation, whereas for multiaxial, the different 
pixels of the detector scans different OPDs. The nature of the recorded signal is 
then different and different biases affect them. 

The Van-Cittert / Zernike theorem: 



The Van-Cittert / Zernike theorem describes the relation between what we 
call the "complex visibility" (i.e. /ii^2 = 71,2(0)) of an object and its brightness 
distribution o{a) on t he plane of the sk y. One can find its demonstration in 
various books, such as iGoodman I (119851 ). 



Theorem: For a non-coherent and almost monochromatic extended source, 
the complex visibility is the normalised Fourier transform (hereafter FT) of 
the brightness distribution of the source. 

By definition, the visibility /i is a complex number, whose modulus is between 
and 1 (see Fig. 9 for an example in the case of a mono-pupil instrument). It 
is by definition: 



fi{li 



o[ u 



no{is) 



(10) 



where ~ represents the FT. By definition, when one measures ||/i|| < 1, the 
object observed is resolved by the instrument. We can see here that interfer- 
ometry does not provide direct access to the image of an object, as single-dish 
telescopes do. It is sensitive directly to the FT of the 



brightness distribution of the object. 



1.6 Practical considerations 



Optical long-baseline interferometry has several intrinsic difficulties: 

(1) the intrinsic complexity to measure the data, 

(2) the apparent complexity for understanding the interferometric measures, 

(3) the very sparse sampling of data, 

(4) the relative bad sensitivity compared to single-dish experiments. 



1.6.1 Measurement complexity 

The ffist problem can be divided into two main problems: 

• The atmospheric effect on ground-based observatories produces very fast 
and varying phase variations. In the case of single-dish telescopes, this ef- 
fect typically reduces the angular r esolution to the one of a telescope of 



the Fried diameter tq (JFried I . Il965l ). Today, the advent of adaptive optics 
allows one to reach much higher angular resolutions at the cost of monitor- 
ing the atmosphere effects very fast and limiting the accessible magnitude. 
In the case of interferometry, long exposure times are not possible due to 
this atmospheric phenomenon. Therefore, just like with adaptive optics, a 
fast feedback loop must operate to correct the fringes motion. This is the 
only way of getting longer exposure times with the science instrument (see 
Fig. 5). 

The use of many subsystems in an interferometer (telescopes, delay lines, 
image sensors, fringe sensors, adaptive optics, focal instruments, etc.) can 
also be an issue. The more complex the interferometer gets, the higher the 
probability of failure. System optimisation is then a key point to manage 
these complex systems. 



1.6.2 Interferometry understanding 

The second problem is about to be solved with the development of general 
astrophysics instruments (AMBER is an example), and the preparation and 
interpretation of observations can be done with easy-to-use tools (developed 



Incident wavefiont 
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with different temperatures 
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Adaptive optics 
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Fig. 5. This figure shows the main effects affecting an interferometric measure- 
ment: the atmospheric OPD being the most dominant by far. It also shows today's 
solutions to correct these effects: adaptive optics (or a tip/tilt correction for small 
apertures) for telescope effects, and fringe tracker and delay lines for interferometer 
effects. 

by ESOilJ, JMMCLiI, and MSCLiI). For the general user, continuous trainings 
and summer schools are intended to help with the data interpretation. 



^ European Southern Observatory 
^ Jean-Marie Mariotti Centre 
^ Michelson Science Centre 



10 



1.6.3 Sparse sampling 

The third problem is also a 2-part problem: 

• The limited number of telescopes (currently 3 on the VLTI, 4 with the sec- 
ond generation, 6 maximum for CHARA) generates a very sparse sampling 
of the (u,v) plane. 

• The atmosphere makes it very difficult to access the object's phase and pre- 
vents a direct measurement of it. Several methods exist to partially retrieve 
the phase: 

■ the phase closure (starting from 3 telescopes), 

■ the differential phase (using a spectrograph), 

■ the phase reference. 

Today, the phase closure is the most widely used since it is easiest to cali- 
brate. 

The strategy to observe and interpret the data is then very different from what 
can be done with imaging facilities (such as single-dish telescopes, at a lower 
spatial resolution). 



1.6.4 Low sensitivity 

The last problem is also about to disappear with the advent of new generation 
interferometers like Keck-I or VLTI, combining long baselines and very large 
apertures. 



2 (u,v) plane properties 



As was seen in the previous section, an interferometer is sensitive to the FT 
of an object's brightness distribution. The question of Fourier plane (also 
called (u,v) plane) sampling is then crucial to know what part of the object's 
information the observer really observed. This section explains how to get a 
sufficiently good (u,v) coverage for an observation using a given interferometer. 



2.1 Super- synthesis 



The idea behind this word is to use Earth rotation to get a larger (u,v) sam- 
pling. Since a stellar interferometer baseline is fixed on the ground, its projec- 
tion on the sky plane changes as the Earth rotates. This baseline projection 
depends only on the hour angle h (i.e. h = LST - R.A.), the baseline coordi- 
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nates (X, Y, Z) and the declination of the object S. A change of coordinates 
hnks the basehne position and its projection on the plane of the sky: 



/.^ 



V 



1 
A 



sin{h) cos{h) 

—sin{6)cos{h) sin{6)cos{h) cos{6) 
cos{6)cos{h) —cos{5)sin{h) sin{6) 



\fx\ 



Y 



:iii 



The (u,v) coordinates for a given object (5 fixed) depend only on a linear 
expression of cos{h) and sin(/i). These coordinates lies on an ellipse, as seen 
in Fig. 6. The time spent on the object will determine the increase in (u,v) 
coverage. However, the relation between the (u,v) filling and this time span 
cannot be intuitively estimated. Therefore, one has to check this for each 
studied case. 



2.2 Number of telescopes 



The relation between the nur nber of tel escopes combined and the number of 
baselines is the following (See iMonnier 1 . l2003i) : 



Nh 



NtiNt 



(12) 



The quantity measured is the amplitude and phase of the visibility function 
at the given baselines. Therefore, the total number of measurable information 



is: 



N: 



tot 



NtiNt 



(13) 



However, in optical long baseline stellar interferometry, things are not that 
simple. When using one spectral bin (for example, with a large band instru- 
ment) the phase information on each baseline is completely lost due to the 
atmosphere random phase noise insertion. However, some phase information 
can be retrieved with the phase closure when using 3 or more telescopes. The 
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Fig. 6. Starting from a snapshot with 2 telescopes (UT1-UT4, top left), and 
increasing both the time sampling (1.5h top-right and l/2h bottom-left) and the 
total integration time (1/2 night top-right and 1 night bottom-left), one can see 
that the (u,v) path is an arc of an ellipse (bottom-right). 

number of phase closures measurable is: 



Nr 



{Nt - 1) {Nt - 2) 



(14) 



^Mi|^ amplitudes and i^X-lM*^ 
interferometry. This gives the actual measurable number of information: 



One normally has access to "*''^^' ''' amplitudes and '^"' '''^""^ ^' phases in 



7V-* = (AT, - 1)^ 



(15) 
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The comparison between eq.l2 and 14 (see Fig. 7) and also between eq.l3 
and 15 suggest that the more telescopes you have, the better you can measure 
information at once. 



1000. 




# Telescopes 

Fig. 7. One gain on 2 sides when adding telescopes to an optical long baseline stellar 
interferometer. The first is the squared increase of available baselines (dashed line). 
The second advantage is the phase closure information (full line) that grows faster 
for a few number of telescopes and tends to the maximum number of measurable 
phases (one per baseline, dashed line) for a high number of telescopes. 



2.3 Spectral coverage 



Spectral coverage helps in filling the (u,v) plane simply because different wave- 
lengths probe different spatial frequencies for a given baseline: 



B 



(16) 



Thus, using a spectroscopic interferometer allows one to scan a line in the 
(u,v) plane with only one snapshot (see Fig. 8, top-right). The more spectral 
range an instrument has, the more (u,v) plane filling it gives. 

Moreover, spectral coverage also allows one to greatly improve the number 
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Fig. 8. Top-left: (u,v) coverage of a single snapshot with a single spectral bin 
instrument (ex: VINCI). Top-right: The same but using spectroscopic capability 
(ex: AMBER 2T). Middle- left: Adding super-synthesis to the last case improves 
(u,v) coverage. Middle-right: The same as before but using 3 telescopes (AMBER 
today). Bottom-left: Using 4 telescopes. Bottom-right: Using 8 telescopes (4 
Ats, thin lines and 4 UTs, thick lines). 



of measures one can access since more p hase info r matio n is accessible via 



differential phase, as explained in detail in iMillour I (120061 ). 



However, such improvement assumes that the object's shape is achromatic 
with regards to wavelength. For example, in the specific case where the object's 
shape is different for each wavelength, then this spectral coverage property 
does not improve (u,v) coverage. 

As a matter of conclusion, combining many telescopes, a long observing time, 
and spectral coverage allows one to get a better (u,v) coverage of an object. 
Since the number of telescopes is very limited and many instruments were 
not able to spectrally analyse the light until recently, the only way was to 
increase the exposure time. Today, both the spectral coverage and the number 
of telescopes is increasing (6 telescopes for CHARA and simultaneous J, H and 
K bands for AMBER are examples), which allows one to begin imaging of the 
targets. 



3 The shape-visibility relation 

Here I add several notations to the previous ones: 



• P=\\u\\, 

• a,l,L are the object's typical dimensions (diameter, FWHM, etc. if any), 

• i? is a flux ratio (if any). 

As seen in previous sections, the measurement of an interferometer is related 
to the FT of the brightness distribution of the object. The goal of this section 
is to show the (complex) visibility description and to have some generic laws 
one can use to perform a simple interpretation. 



3.1 Generic properties 



First of all, one always has to have in mind that interferometry deals with 
FTs. We recall here the generic properties of FTs for continuous functions: 

• linearity (addition): FT[f + g] = FT[f] + FT[g], 

• translation (shift): FT[f{x -xo,y- yo)] = FT[f]{u,v) x e2»'^'("^o+''w)^ 

• similarity (zoom and shrink): FT[f{ax, Py)] = ^FT[f]{^,^), 

• convolution ("blurring"): FT[f ^ g] = FT[f] x FT[g], 

• oo limit ("small" details): FT[f] ^^ 0, 
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• limit ("large" details): FT[f] ^ 1. 

The last two points lead to a specificity of optical long-baseline stellar interfer- 
ometry: the best constraints for a model will be given by a measured visibility 
of 0.5. This translates into: 

• For short baselines (i.e. small spatial frequencies), all the possible shapes 
have degenerate visibilities. The visibihty modulus has a squ ared depen- 



dency with base length (see Fig. 10 and the demonstration in iLachaume 



20031 ) and its phase has a linear dependency. 

• For long baselines, all continuous models (i.e., the ones not including S 
functions) have visibilities towards 0. For a given accuracy on the measure- 
ments, therefore, no discrimination between different models can be done 
as the object is "over- resolved" . 

Therefore, to maximise the efficiency of an observation, one has to know the 
approximate shape and size of an object in advance (1^* guess). The knowledge 
of an object must then be as high as possible to maximise the chances of success 
of an interferometric observation. Once the observations have been made, and 
for a given model, there are many ways to fit the visibilities without changing 
the model itself: 

• change the distance (using similarity), 

• change the orientation in the plane of sky (using linearity and coordinate 
transformations) , 

• shrink or expand the model on one axis (using similarity), or 

• blur the model (using convolution). 

The two first points are useful for model fitting as one model image FT can 
be scaled and rotated to fit the observed data. 



3.2 Examples 



Here 1 present typical examples where the visibility function can be computed 
easily with an analytical formula. They are the basis of all further analysis, 
allowing one to get a first idea of the object's shape without doing a complete 
physical modelling. 



3.2.0.1 Point source: double or multiple star A centered point source 
is described in the Fourier plane by a constant. For many stars, even if the re- 
solving power of interferometry is much higher than for single-dish telescopes, 
a good approximation of their shape can be made, assuming they are unre- 
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solved. The shift and add properties of the FT give the visibihty expression 
for a multiple star: 



V{lt\ 



ELihcos{^) 



^N 



Efc=i 4 



(17) 





U(m) 



U(m) 



Fig. 9. Example (taken from iBalega et al. 1 . 12006| ) of a diffraction-limited image 
of a binary star (HIP 4849, observed at the Special Astronomical Observatory, 
Zelentchouk, left). The complex visibility of this image (middle and right) gives 
typical waves until the cut-off frequency of 6 m (features at larger frequencies are 
produced by the image reconstruction algorithm used here). 



3.2.0.2 Gaussian disk: extended envelope In a first approximation, a 
Gaussian brightness distribution can describe several types of envelopes, such 
as an opaque wind surface around Wolf-Rayet stars, gaseous or dusty disks 
around young stars, etc. A Gaussian brightness distribution has a Gaussian 
visibility function: 



Vip) 



(Trap)_ 



3.2.0.3 Uniform disk: stellar surface In a first attempt, a uniform disk 
brightness distribution is well suited to describe a stellar surface, since, in gen- 
eral, the stars look like a sharp-edged disk, corresponding to the photosphere 
surface. The visibility function of a uniform disk is a 1'^* order Bessel function 
divided by Trap: 



Vip) 



Ji(7r ap) 
nap 



(19) 




Fig. 10. Whichever model is used (binary star in solid line, Gaussian disk in dotted 
hne, uniform disk in dashed hne, or a composite object in dash-dotted line), the 
smaU basehnes do not allow one to distinguish it's shape. 



3.2.0.4 Ring and Hankel function: any circular object A ring can 
describe a thin shell around a star. Its expression depends on the 0*^ order 
Bessel function: 



V{p) = Jo(vrap) 



(20) 



It is also very interesting to use the addition property of the FT to produce 
the visibility function of any circular object with a given profile /(r): 



V(p) = 2tx \ I(r)3o(2nrp)rdr 
Jo 



(21) 



3.2.0.5 Other objects For other objects, however, no simple tool except 
FT can be used. The best thing is to produce an image for the object in 
each wavelength of interest and then produce its FT and compare it with the 
visibilities. 

In Table 2, one can see a summary of the different possibilities offered to test 
simple models. 
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Table 2 
Summary table for simple shapes and their corresponding complex visibility. 



Shape 



Brightness distribution 



Visibility 



Point source 



d{^ 



1 



Binary star 



A[5{T}+R6{^-J^)] 



l+i?2+2_Rcos 



T+W 






Gauss 



41n(2a) 



X e 



-4 In 2 



7^ 



(Trap) 
e 4 In 2 



Uniform disk 



otherwise 



Jl(7rap) 
Trap 



Ring 



Tra V 2/ 



Jo (Trap) 



Any circular object 



/(r) 



27rL I{r)3o{27rrp)rdr 



Pixel (image brick) 



4- if X < / and y < L 



IL 

otherwise 



sin(7ra;i) sin{7rj/L) 
TT^xylL 



4 Conclusions 



I have reported the main things to know about astronomical optical interfer- 
ometry: a short history, the main principles, and the main properties of the 
measurements one gets from optical long baseline interferometry. Using this 
information, and by doing the practical exercises associated with this article, 
a future observer should be able to prepare observations knowing the real 
limitations but also the true possibilities of today's interferometers. 
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